Fluctuation Theorem for the flashing ratchet model of molecular motors 
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Molecular motors convert chemical energy derived from the hydrolysis of ATP into mechanical 
energy. A well-studied model of a molecular motor is the flashing ratchet model. We show that 
this model exhibits a fluctuation relation known as the Gallavotti- Cohen symmetry. Our study 
highlights the fact that the symmetry is present only if the chemical and mechanical degrees of 
freedom are both included in the description. 
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Molecular motors are subject to intense study both 
from biological and technological point of view [l|, 
These remarkable nanomachines are enzymes capable of 
converting chemical energy derived from ATP hydrolysis 
into mechanical work. They typically operate far from 
equilibrium, in a regime where the usual thermodynami- 
cal laws do not apply. Generically such motors are mod- 
elled either in terms of continuous flashing ratchets [1, 01 
or by a master equation on a discrete space [1, Q . Re- 
cently, a general organizing principle for non-equilibrium 
systems has emerged which is known under the name of 
fluctuation relations These relations, which hold 

for non-equilibrium steady states, can be seen as macro- 
scopic consequences of generalized detailed balance con- 
ditions, which themselves arise due to the invariance un- 
der time reversal of the dynamics at the microscopic scale 
i- 

An interesting ground to a pply these concepts is the 
field of molecular motors [lo, [ill, [11, III [11 [11, III. 
The fluctuation relations impose thermodynamic con- 
straints on the operation of these machines, particularly 
in regimes arbitrary far from equilibrium. Near equilib- 
rium, they lead to Einstein and Onsager relations. For 
non-equilibrium steady states, they can be used to quan- 
tify deviations from Einstein and Onsager relations as we 
have shown in Refs. (isl. [l3|. 

In this paper, we investigate fluctuation relations for 
continuous ratchet models. We flrst study a purely me- 
chanical ratchet (model I) , which applies to the transloca- 
tion of a polymer through a pore [l7| . We then consider 
a flashing ratchet (model II), which applies to molecu- 
lar motors [3]. Using a method inspired by Refs. [3, [11], 
we show that the Gallavotti-Cohen symmetry is always 
present in model I, but we emphasize that in model II 
the symmetry exists only if the chemical and mechanical 
degrees of freedom of the motor are both included in the 
description. 

Let us flrst consider a random walker in a periodic 
potential subject to an external force F (model I) [^.[igj. 
The corresponding Fokker-Planck equation is 



dt ^ dx 



dP 
dx 



U'jx) ~ 
kuT 



where U{x) is a periodic potential U{x + a) = U{x) and a 
is the period. This equation describes the stochastic dy- 
namics of a particle in the effective potential Ueff{x) — 
U{x) — Fx. B y s olving Eq. ([T]) with periodic boundary 
conditions [TtI. [T9j. it can be readily proven that the sys- 
tem reaches a stationary state with a uniform current J 
in the long time limit. This current is non- vanishing if a 
non zero force is applied. When F = 0, there is no tilt 
in the potential, J — and the stationary probability is 
given by the equilibrium Boltzmann-Gibbs factor. 

We call x{t) the position of the ratchet at time t know- 
ing that the ratchet was located at x{0) =0 at time t = 0, 
which wc decompose as x — (n + (^)a where n is an inte- 
ger and < C < 1. The stationary current J is related 



to the average position x(t) hy J = limj 



{x{t)) 



i.e. J 



is the mean speed of the ratchet in the long time limit. 
More generally we are interested in the higher cumulants 
of x{t) when i — > oo. It is useful to deflne the generating 
function 

Fa(C, t)^Y. c'^P W< + ")) ^(("^ + Oa, t) . (2) 

n 

The time evolution of this generating function F\ is ob- 
tained by summing over Eq. ([ij. This leads to the fol- 
lowing equation: 



at 



C{X)F^iU)^ 



(3) 



where the deformed differential operator C{X) acts on a 
periodic function $(C,i) of period 1 as follows: 



(9$ 

(4) 

where t/g^y = adxUef f /ksT and the left hand side of 
Eq. (|4]) is proportional to the inverse of the characteristic 
time T = a? I Dq. A similar procedure exists in solid 
state physics, where periodic functions are expanded in 
eigenfunctions of Bloch form, which are eigenfunctions of 
an operator similar to £(A) [l7| . 

The operator C{X) has the following fundamental con- 
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FIG. 1: Sketch of the sawtooth potential U{x). The poten- 
tial has period a, aa is the distance from a minimum to the 
next maximum on the right, and Uo is the maximum of the 
potential. 



jugation property: 

^Ui.)/ksT^(^^^ (^^-Ui.)/ksT^'^ ^ (_ J _ X) (5) 

with / — FalksT the normalized force. This property 
implies that operators £(A) and (— / — A) are adjoint 
to each other, and thus have the same spectrum. If we 
call G(A) the largest eigenvalue of >C(A), we obtain from 
Eq.[n]that 0(A) satisfies the Gallavotti-Cohen symmetry: 

e(A) = e(-/-A). (6) 

In fact, this symmetry holds for all eigenvalues. For the 
special case / = 0, the conjugation relation ([5]) reduces 
to the detailed balance property [isj . Finally, it is impor- 
tant to note that 0(A) is the generating function for the 
cumulants of x{t). 

We have calculated numerically the function 0(A) for 
the case of the sawtooth potential shown in F ig. [H with 
a barrier height Uq of order of several fc^T This 
function was obtained by first discretizing the operator 
£(A) and then calculating its largest eigenvalue using the 
Ritz variational method. This method does not require to 
find a basis specific to the chosen potential, in contrast to 
what was done in Ref. 20] for the cosine potential. Our 
numerical method can handle any shape of the potential. 

The form of 0(/ry) with -q = A// is shown in Fig. [2] 
for different values of the normalized force /. The 
symmetry of all the curves with respect to rj = 1/2 
corresponds to the symmetry of Eq. At weak force, 
Q{frj) has a parabolic shape associated with gaussian 
fluctuations, whereas at higher forces a flattening occurs 
associated with non-gaussian fluctuations [l^, 13, HO]- 
By numerically taking derivatives of 0(A) with respect 
to A near A = 0, we recover the velocity obtained by 
directly solving Eq. ([1]) [l3,[ll. 

We now come to the derivation of the Gallavotti-Cohen 
symmetry for the flashing ratchet model (model II). In 
this model [1, [2l|, (2^ , the motor has two internal states 
i — 1,2, which are described by two time- independent 
potentials Ui{x). We assume that these potentials are 
periodic with a common period a. The probability den- 
sity for the motor to be at position x at time t and in 
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FIG. 2: Normalized eigenvalue TQ{frf) (with r = a'^ /Do) as 
function of rj for different values of the normalized force /; 
from top to bottom, / = 5, / = 10 and / = 20. The parame- 
ters of the potential are a = 0.7, Uo/kBT — 5. The symmetry 
of all the curves with respect to r; = —1/2 is Gallavotti-Cohen 
symmetry expected for model I. 




FIG. 3: Normalized eigenvalue TQ{frj) as function of r; for a 
normalized force f — 5 (top two curves) and / = 10 (bottom 
two curves) for the fiasing ratchet (model II). The solid curves 
correspond to the case where the switching rates satisfy de- 
tailed balance which leads to the Gallavotti-Cohen symmetry. 
The curves with circles (/ — 5) and crosses (/ = 10) corre- 
spond to cases where detailed balance is broken with constant 
switching rates uji{x) — uj2{x) — 10r~^ and with the same po- 
tentials. The lack of symmetry in these curves with respect 
to = —1/2 is apparent specially near 77 = —1. 

state i is Pi {x,t). The dynamics of the model is described 
by 

^ + ^ = ~U;,{X)P,+LU2{X)P2 

r)P r) T 

^ + ^ ^ cj,{x)P, ~ u;2{x)P2 , (7) 
where lui{x) and lu2{x) are space dependent transition 
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rates, and the local currents Ji are defined by: 



J,; = -De, 



dP,, 



1 



dx knT V dx 



dU,, 



(8) 



with Dq the diffusion coefhcient of the motor and F a 
non-conservative force acting on the motor. The tran- 
sition rates can be modeled using standard kinetics for 
the different chemical pathways between the two states 
of the motor [2l| : 

LU2ix) - -f V(x)]e(^^(")-^")/'^'«^, (9) 

where A/i = Afl/ksT is the normalized chemical po- 
tential and A/i the chemical potential associated with 
ATP hydrolysis. Terms proportional to uj{x) are associ- 
ated with thermal transitions, while terms proportional 
to ^"(2;) correspond to transitions induced by ATP hy- 
drolysis. One could easily introduce more chemical path- 
ways than the ones considered here ^21j but this extension 
is not essential for the present argument. Note that the 
way the force enters the rates is unambiguous in such a 
continuous model 5, fij. 

Note that Eq. ([7]) can be rewritten as a matrix C of 
operators: 



d_fp, 

dt [ P2 



= c 



Pi 
P2 



Pi 
P2 



(10) 



where the action of the operator d on a function $(a;, i) 
is given by 



' dx"^ 



Do — 



d fU'-F_ 



dx \ ksT 



(11) 



When F 
and 



and A/i = 0, the system is at equilibrium 



UJl{x) 



exp 



U2-U1 

knT 



(12) 



In this case, the stationary solution of the system ((7|) is 
the Boltzmann distribution for Pi and P2, the currents 
Ji and J2 vanish and there is no global displacement of 
the motor. If both F and A/j, do not vanish, then the 
system is out of equilibrium and non-vanishing currents 
can appear. 

If the switching between the two potentials occurs only 
by thermal transitions, i.e when A/z = 0, the rates satisfy 
the detailed balance condition of Eq. (fT2|) . even in the 
presence of a non-zero force F. The Gallavotti-Cohen 
symmetry follows by considering a 2x2 diagonal matrix 
of operators Ci{X) of the form The symmetry is in- 
deed present as shown in the solid curves of Fig. [3] In 
the general case however, where the normalized force / 
and chemical potential A/i are both non-zero, the rela- 
tion is no more satisfied and the Gallavotti-Cohen 
relation ^ is not valid. This is shown in the curves with 
symbols in Fig. [3] where for simplicity we took constant 




FIG. 4: For model II, the normalized eigenvalue 
rA(/r/, — A/i/2) is shown as function of rj. The dashed curve 
corresponds to / = 5 and A/i = 0, the solid curve corresponds 
to / = 5 and A// = 10, the symbols correspond to f — 2 and 
A/i = 10. The symmetry is recovered in all cases in this de- 
scription which includes both the mechanical and chemical 
degrees of freedom. 



switching rates wi = W2 = 10t~^. For all the curves 
of this figure, we took a sawtooth potential Ui with the 
same parameters as in Fig. [2l and a potential U2 con- 
stant in space. The breaking of the symmetry of Eq. [5] 
can be interpreted as a result of the existence of internal 
degrees of freedom, similarly to the violations discussed 
in Ref. [H. 

To establish a fluctuation relation for the flashing 
ratchet model, one must consider both the mechanical 
and chemical currents present [l^ [l^. Let us intro- 
duce the probability density Pi{x,q;t) associated with 
the probability that at time t the ratchet is in the in- 
ternal state i, at position x and that q chemical units of 
ATP have been consumed. The evolution equations for 
this probability density is obtained by modifying equa- 
tions ([7]) after taking into account the dynamics of the 
discrete variable q. We have 



dPijx^qA) 
dt 

dP2{x,q,t) 
dt 



= {Ci-LUi{x))Pi{x,q,t) (13) 

+ uj^\x)P2{x, q + l,t)+ LJ°{x)P2{x, q, t) 

= iC2-i02{x))P2ix,q,t) (14) 

+ u;°{x)Pi{x,q,t) + u;l{x)Piix,q - l,t). 



We use a notation similar to that of Ref. [l3|, where w-(a;) 
denotes the transition rate at position x from the internal 
state i with I = —1, 0, 1 ATP molecules consumed. This 



-fx)/kBT^ ^1 



leads to = uje^u^-fx)/kBT ^ ^0 ^ ^^(u^ 

i^i{x) + culix) and uj2{x) — uj2{x) + uj2^{x). 

As above we introduce two generating functions -F'i.a,7 
and F2,A,7, depending on two parameters A and 7 which 
are conjugate variables to the position of the ratchet and 
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to the ATP counter q. We have for i = 1,2, 
F.,x,^iC = E E e^^'^+"^^.(«(C + n), q; t). (15) 

g n 

The evolution equation for these generating functions is 
obtained from Eq. p5|) as 

with the operator £(A,7) decomposed as 

£(A,7)-2?(A)+AA(7), (17) 

with 2?(A) the diagonal matrix diag(£i(A) — lui, C2{X) — 
UJ2), and 

Consider now the diagonal matrix Q defined by 
diag(e~'^^/'^^^, e~^'^/^'^'^). By direct calculation, one can 
check that Q-^Af{j)Q = Af^-A^i - 7) . From Eq. ©, 
one obtains Q~^'D{'-f)Q — 2?^ (— A/i — 7). By combining 
these two equations, we conclude that 

g-i/:(A,7)Q = /:U-/-A,-AM-7), (w) 

which leads to the Gallavotti-Cohen symmetry: 

A(A,7) = A(-/-A,-A/i-7), (20) 

where A(A,7) is the largest eigenvalue of >C(A,7). If we 
consider only the mechanical displacement of the ratchet, 
the relevant eigenvalue 6(A) is given by 6(A) = A(A, 0), 
which clearly does not satisfy the fluctuation relation 
as shown in Fig. [S] In Fig. [4l we have computed 
A(/77, — A/z/2) for the same potentials and with rates 
i^l{x) of the form given above with lli{x) ~ hT~^ and 



4>{x) = lOr ^. We have verified that in all cases the 
symmetry of Eq. holds. 

In this paper, we have shown that the large devia- 
tion function of the mechano-chemical currents obeys the 
Gallavotti-Cohen relation. Another related but different 
symmetry relation for the entropy production exists un- 
der more general conditions [tI [lOlTlsl [Tsl. [20|. We have 
shown here that the symmetry for the currents is valid for 
the flashing ratchet model when internal degrees of free- 
dom are taken into account. This raises a fundamental 
question concerning the validity of fluctuations relations 
and their applicability to other types of ratchet models 
0, 01 ■ Other mechanisms exist which are known to pro- 
duce deviations from fluctuations relations [23| . and it 
would be interesting to investigate whether fluctuations 
relations can always be restored by a suitable modifica- 
tion of the dynamics. 

On the experimental side, it would be very interesting 
to investigate fluctuations relations for molecular mo- 
tors using single molecule experiments, in a way simi- 
lar to what was achieved in colloidal beads or biopoly- 
mers experiments 0. Using fluorescently labeled ATP 
molecules, recent experiments with myosin 5a and with 
the Fq — Fi rotary motor, aim at simultaneous recording 
of the turnover of single fluorescent ATP molecules and 
the resulting mechanical steps of the molecular motor 
[23 |. These exciting results indicate that a simultaneous 
measurement of the values of the mechanical and chemi- 
cal variable of the motor is achievable, and therefore from 
the statistics of such measurements it is possible to con- 
struct P{x, q, t). With enough statistics of such data, one 
could thus in principle verify Eq. 1201 Such a verification 
would confirm that the Gallavotti-Cohen symmetry is a 
thermodynamic constraint that plays an essential role in 
the mechano-chemical coupling of molecular motors. 
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